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ON THE MOTION OF THE CENTRE OF QRA VITY OF 
THE EARTH AND MOON. 



That the motion of the centre of gravity of the earth and moon is sen- 
sibly the same as if the masses of these two bodies were concentrated at this 
centre has been confidently asserted over and over again. However a little 
scepticism on the matter may not be altogether ill-advised. Were this as- 
sertion true it would follow that, setting aside the action of the planets, we 
should get the sensibly exact mean angular motion of this centre about tha 
sun, by first deriving the mean distance a' from the elliptic value of the 
radius vector 

r' = a'[l-fje' 2 -f periodic terms], 

and then n' from the equation n' = -J-y^, M denoting the sum of the 

masses of the sun, earth and moon. 

Let us see whether this value is sensibly exact under the conditions we 
suppose. 

Assume that the masses of the sun, earth and moon are denoted by wtj, 
m 2 and m 3 , and their rectangular coordinates severally by $ lf tj 1} £ x ; £ 2 , 
Vi> C2 > ^3> Vs> £s- And let the rectangular coordinates of the moon rela- 
tive to the earth be denoted by x, y and 2; those of the : 
sun relative to the centre of gravity of the earth and 
moon by x', y' and z' ; and those of the centre of grav- 
ity of the three bodies by X, Y and Z. Then from ! 
an attentive consideration of the subjoined figure,where ' 
S, E and M denote the positions of the sun, earth and 
moon, C the centre of gravity of the last two bodies, 
and C the centre of gravity of all three, it will be seen 
that, if we put 
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u = m s „r _ m 2 +m 3 

m 2 +m s ' ' m 1 +m 2 +m 3 ' 

we shall have 

$ 1 = f z'x'+X, 

f 2 = Gm'-ik— /«+x, 

with two groups, of three equations each, for the -q and £, obtained from 
these by writing, in the second members, for x and X, y and Y } and again 
2 and Z. 

If we differentiate the equations just written, then square and add the 
results after having multiplied them severally by m t , m 2 and m 3 we shall 
get 

m^l+m^l+m^l = m 1 ft'dx ,2 +m 2 fidx' i +MdX i . 

From this equation it is evident that, if Q denote the potential function, 
the differential equations, determining the variables x, y, z, x', y', z', are 



d 2 x dQ 
m ^W = d^> 


,dV dQ 
m ^-W=dlc» 


d*y dii 


,d\' dQ 
m ^-d¥ = dy-» 


d*z dii 


.dh' dQ 
ir dt 2 dz' 



Hence it may be gathered that the disturbing function for the motion of 
the sun relative to the centre of gravity of the earth and moon differs from 
the corresponding function for the motion of the moon relative to the earth 
only by a constant factor which depends on the masses. 

The expression for Q is 

O = m l m 2 I m l m 3 I m 2 m 3 
4L2 + ^1.8 + ^2,3 ' 

where the J's are given by the equations 

^ >2 = (x'+fixf+(y'+fiyf+(z'+fizY, 

4l 3 =x*+y*+z\ 
Let us put 
r» = ar'+y+z 2 , r n = x' 2 +y n +z n , rr'S = xx'+yy'+zz'. 
Then 

A\ 2 =r' 2 +2//rr'/Sf+^ 2 r 2 , 

A \,z = »"' 2 — 2(1— fi)rr , S+(l— fift*. 
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r 1 1 

Since the ratio — T is only about — — , and « about — , it is convenient to 
r' J 400 r 80' 

expand, in ii, the reciprocals of J 12 and J 1>3 in infinite series proceeding 
according to ascending powers of — r . This, in both cases, evidently de- 
pends on the developement of 

(1— 2ax+o?)->i 
in powers of a. By the Theorem of Lagrange, in solving the equation 
y — al\y) = x with respect to y, we get 

, rv x , « 2 d.F(xY , , a" d"-i.F(x) n , 

whence 

dy _. dJ^.a* d?JJ\ocf , a n d n .F(x) H _ L 

Let us suppose that we have here F(y)= Ky 2 — 1); the equation, on which 
y depends, becomes then y — ^(y 2 — 1) = x, and the resolution of this 
quadratic in y gives 1 — ay = -j/(l — 2aa; -f- a 2 ), and by differentiation 

^ = (1— 2ox+a 2 p a \ Consequently 

+a [jA 1'2AJ 

3 r6-5.4^ 3 4.3.2 -] 

+a |_2A6 I" 2.4.6 * J 

^ r 8.7.6.5 4 __4 6.5.4.3^,4.3 4.3.2.1 H 

^ L2.4.6.8 1' 2.4.6.8 T 1.2* 2.4.6.8 J 

+ 

The law of the numerical coefficients in this series is so plain that we can 
set down as many terms as we have occasion for. 

In making the application to the reciprocals of d 1 2 and J 13 we must 

put, in the first case, a — — fi-„ in the second, a = (1 — //)-, and in both, 

x = S. 

We obtain as the potential function proper for the relative motion of the 
moon about the earth 

1 n m, +m<> . fl 1 i 1 1 ~1 
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__m 2 +ro 3 
r 

+m 1 {[(l-^)-i+^]l 

+ >■ 

To get the similar function for the relative motion of the sun about the 
centre of gravity of the earth and moon, it is necessary to muliply the pre- 
ceding expression by 

m 2t* _ «h+^2+»*3 /( (l_ ;; ). 



The term of the potential function for the moon, factored by r 3 -*-/ 4 , 
gives rise to inequalities in the lunar coordinates factored by a-^-a'. As 
this term has 1 — 2/x as a factor, we see the correctness of the rule which 
directs to multiply this class of inequalities by 1 — 2fi in order to include 
the effect of the disturbance of the relative motion of the sun about the earth 
by the lunar mass. 

In treating the motion of the sun about the centre of gravity of the earth 
and moon it will suffice to take two terms of the preceding expression and 
put 

Let the longitudes of the sun and moon be denoted respectively by X' and 
X, and neglect the latitudes; then 

^0= ^ +i ifKl-^[3cos2(W)+l]. 

The differential equations, determining r' and X', are 

% -r'% +^V*l-rt£[8 cos 2r+l] = 0, 

U r 'w) -f^'W-^sin2r =0, 

where, it will be noticed, we have put r = a, and, after differentiation, in 
the final small terms, r'—a', M-±-a' 3 = n' 2 and X — X = t the mean angular 
distance of the moon from the sun. The integration of the second eq'n gives 
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a W 3 i 

dt 



dX aJn' 3 n n n .a 2 „ 
r n 4n — n' v 'a* 



a being the arbitrary constant. We can now eliminate — ,— from the first 

equation and we get 

dV , M ain' 2 . , ,, , ,, -.a 2 r3n— n' o . i ~1 n 
-ST+-72 & 7r +|w'W//(l— //)-^ -cos2r+l = 0. 

Let us suppose that this equation is satisfied by 
r' = a -^a'a^ cos 2r, 
a a being a coefiicient to be determined. Substituting this value of r' in the 
differential equation we get the two equations of condition 
M 



i% 



•n n a +tn n a't41—ftF n = 0, 



(4n 2 — 8nn'+3n' 2 )a,— |n'2^=^-«(l— //)— = 0. 
1 n — n' a' 2 



Whence may be derived 



2 



^ = 1+^(1-^, 

__ 3 m 2 (3— m) n— \ a * 

ai ~ 4 (1— m)(4— 8m+3m 2 ) A W a' 2 ' 
where, as is usually done in the lunar theory, we have put n'-i-n—m. The 

value of r r , thus obtained, being substituted in the expression for -j-, we get 

dA' , , , m 4 — 2wi+m 2 ,., .a 2 

_ = «'— fn/ — - — J-g-j fi{l—fjt)-r= cos 2r. 

eft 1 — ro4 — 8m+3m 2 ^ v ■'a 2 

Integrating 

v / i u 3 m 2 4 — 2m4- m 2 /.. N a 2 . 

A' = e'+n'i — -r — -^ - — - — \-~— 2 fjt(l— fA-TzSm 2r. 

4 (1 — m) 2 4 — 8m+3m 2 ^ v ^a 2 

The numerical values of the constant quantities, which enter into these 
formulae, are 

m = °- 0748 ' ^™ > ^= - 002587 ' 

n' = 1295977".4. 
They give us 

r' = a'[1.00000 00200+0.00000 00003 cos 2r], 
X' = e'+n't— 0".0001 sin 2r. 
The periodic terms of these equations are too small for consideration, but 
the constant term of r'-t-a' may be noticed. If we should obtain the value 
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of a' from measured values of r', on the assumption that the value of the 
constant term is unity, it would be too large by the 0.00000 002 part. And 

this value substituted in the equation n' = -y-73, would give n' too small 

by the 0.00000 003 part, or ri would be too small by 0".03895; or the er- 
ror in the mean longitude of the sun would amount to nearly 4" in a century, 
a quantity which could not, in the present state of astronomy, be neglected. 
However, it is only fair to state that astronomers proceed in a way the re- 
verse of this; that is they observe n' and thence deduce a', and in this case 
the term 0.00000 002 is without significance, since the logarithms of the 
radii vectores in the ephemerides are usually given to 7 decimals only. 



A CASE OF SYMBOLIC VS. OPERATIVE EXPANSION. 



BY A. S. HATHAWAY, CORNELL, UNIV., ITHACA, NEW YOEK. 

(d d \ m 

an h«2 j (" •••) > tne 

general term of which is A m a\^ . . f -= — J ( -? — J ..,where.4 m = f t * , 

r-\-s-\- . . . = m; and by {D) m , the operation D or f a^ — \-a 2 -= \- ... J 

repeated m times, whose general term is 

the extra parenthasis here, and in what follows, inclosing a symbol which 
combines operatively. Then (D) 2 differs from D 2 in the production of an 

extra term a< -, — by each portion (D)a.i— — fromf a< -j— J a* -=— of the 

complete operation (D) D or (D) 2 ; so that the complete (D) 2 = D 2 + D. 
And in general, (D) D m ~ l differs from D m by an extra term 

(r+s+...)A^ 1 ^ 2 ...(^J{^'... 

derived from each portion 

{D)A m _ 1 a\a\...(~ — \ (- — )... of the complete operation 

(D)D"-i ; so that, since r+s+. . .=m— 1, (D)D m ~ x =D m +(m—l)D n '^ 
.: D m =[(D-m+l)]-D m - 1 =[(-D-m+l)][(-D~m+2)]..[(I>-l)]-Dor(Z)) , " / . (1) 



